Mixmaster Behavior in Inhomogeneous Cosmological Spacetimes 
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Numerical investigation of a class of inhomogeneous cosmological spacetimes shows evidence 
that at a generic point in space the evolution toward the initial singularity is asymptotically that 
of a spatially homogeneous spacetime with Mixmaster behavior. This supports a long-standing 
conjecture due to Belinskii et al. on the nature of the generic singularity in Einstein's equations. 
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If one assumes that our expanding universe can be described by a spatially homogeneous and isotropic solution to 
Einstein's equations, one can "run the expansion backward" to a hotter, denser universe in the past. Such an analysis 
leads to an understanding of the cosmic microwave background and primordial light element abundances. Running a 
further finite time into the past yields the big bang — a singularity characterized by infinite density, temperature and 
' gravitational tidal force. While this standard cosmological model accounts for observed features of the universe, the 
' reliability of its predictions about features that have not been observed depends on the stability of those predictions 
\ when the conditions of homogeneity and isotropy are relaxed. 

. Einstein's equations allow for a rich variety of cosmological spacetimes, by which we mean solutions that are 
' deterministic (contain a compact Cauchy surface) and have a physically reasonable stress energy tensor (one that 
, satisfies the strong energy condition). Powerful theorems state that, generically, such spacetimes have an initial 
■ singularity. But the theorems do not describe the nature of the singularity. In the approach to the initial singularity 
in the standard cosmological model, the kinetic energy of the isotropic collapse (proportional to the square of the 
ly-j Hubble parameter) dominates the spatial curvature. A similar type of approach to the singularity is found in the 
iy~j Kasner spacetimes Q. These vacuum solutions are anisotropic, spatially homogeneous and spatially flat (type I 
in the Bianchi classification of homogeneous spaces). Since they are spatially flat, the spatial curvature terms are 
absent from the evolution equations for these models, and the kinetic energy of the anisotropic collapse drives the 
approach to the singularity. A cosmological spacetime is said to have an asymptotically velocity term dominated 
(AVTD) singularity if the evolution toward the singularity at each spatial point approaches that of one of the Kasner 
spacetimes or that of a nonvacuum Bianchi I spacetime with fixed Kasner exponents 1^-^ . Another possible behavior 
^ ■ near the singularity is exemplified by the spatially homogeneous Mixmaster spacetimes, in which the dynamics of 
^ ' the collapse to the singularity is approximated by an infinite sequence of Kasner spacetimes with a deterministically 
5-H 1 chaotic transition from one Kasner to the next . In Bianchi type IX the transitions are caused by "bounces" off a 
potential provided by the spatial scalar curvature |^,|^ . In Bianchi type VIq with magnetic field there are bounces off 
^ 1 a potential provided by the magnetic field in addition to the curvature bounces ^ . 

' While much progress has been made in understanding the homogeneous case [^,0, the behavior of spatially 
inhomogeneous solutions to Einstein's equations near an initial singularity is largely unknown. In a number of 
very limited classes of solutions (various classes of spacetimes with 2-torus spatial symmetry and polarized vacuum 
spacetimes with U(l) symmetry) there is strong evidence for AVTD behavior H,|ll|-|l3|. However, few expect AVTD 
behavior to occur generally in cosmological spacetimes. Rather the conjecture has been that generically there is 
Mixmaster behavior, in which the evolution toward the singularity at a generic spatial point approaches that of one 
of the homogeneous Mixmaster spacetimes |6|Jl^ . In a spatially inhomogeneous AVTD spacetime the evolution at 
different spatial points will approach that of different Kasner solutions. In a spatially inhomogeneous spacetime that 
has Mixmaster behavior the evolution at different spatial points will approach that of different Mixmaster solutions. 
Although these two possibilities are mutually exclusive, in both cases the presence of the inhomogeneity ceases 
to govern the dynamics asymptotically toward the singularity. This is a drastic assumption! The space remains 
inhomogeneous at all times, yet the effect of inhomogeneities on the evolution becomes negligible. 

Until now there has been no evidence that Mixmaster behavior occurs in inhomogenous spacetimes. We have 
obtained such evidence, and discuss it in this letter. We have focused on a class of cosmological spacetimes that is 
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an inhomogeneous generalization of Bianchi type VIq with magnetic field. Numerical study of the evolution toward 
the singularity for a representative sample of initial data shows that a regime consistent with Mixmaster behavior 
is reached. It is impossible to follow the evolution all the way to the singularity numerically. But analysis of the 
evolution equations under the conditions that exist late in the numerical evolution show that the regime will continue. 

The spacetime manifold on which these solutions are defined is E"^ x R, where S'^ is a solv-twisted 2-torus bundle 
over the circle [|l5| . While this manifold docs not admit a global two-torus action, it does admit the local group action 
corresponding to Bianchi VIo, which contains a local two-torus action as a subgroup. Spacetimes in the class we 
studied have this local spatial two-torus symmetry and their metrics can be written in the 3-torus Gowdy form 
with appropriate nonperiodic boundary conditions on some of the metric coefficients. In particular, we can write the 
metrics for this class of spacetimes as 
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The metric functions A and fi are periodic in 9. P{6 + 27r, r) = P{9, t) + 27ra and Q{9 + 2ti, t) — Q{9, T)e^^'^". Here a 
is a constant determined by the manifold twist. When A and /i are independent of 9 and P{9, r) = P{9o, t) + a{9 — 9o) 
and Q{9,t) — Q{9o,T)e~°''^^~^°'^ the spacetime is locally homogeneous Bianchi VIq. As in Bianchi VIo with magnetic 
field, we take the Maxwell tensor iohe F = B dx A dy. It then follows from the Einstein-Maxwell field equations that 
B is necessarily constant in space and time. The nondynamical function fi is nonzero only if the electromagnetic field 
is nonzero. The time coordinate r has been defined (without loss of generality) so that the singularity is at r = oo. 

The evolution equations for these solutions of the Einstein-Maxwell field equations can be derived from the Hamil- 
tonian density H — Hi + H2 + H3, 
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H3 = A7rxeS^+-^/^B\ 
In addition, the following constraint equations must be satisfied. 
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npP' + ttqQ' + TTxX' ^0, tta 
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where TTp, ttq and Tr\ are the momenta conjugate to P, Q and A. 

This three hamiltonian form is useful for the following reason. If any one of the three subhamiltonians is taken by 
itself, and the other two ignored, the system is exactly solvable. Thus, after approximating the continuous system by 
a discrete one, Suzuki's decomposition of exponential operators pTf can be used to decompose the evolution operator, 
with each piece exact. The computer code we use for numerical evolution is based on this decomposition. It is an 
adaptation of a code used for the Gowdy spacetimes [ pH , and uses the fourth order decomposition and also fourth 
order accurate representation of the spatial derivatives. 

One can understand the structure of the numerically observed metric evolution for these spacetimes in terms of 
the evolving relative dominance of the three subhamiltonians. Let us fix a spatial point p — {9q, xo,yo). Once the 
Mixmaster regime is reached, then for most of the evolution towards the singularity (r — > 00) at that point, we observe 
numerically that Hi dominates H2 and H3, and the metric evolution is essentially that of some Kasner spacetime. 
Time intervals during which this happens are called Kasner epochs. For intermittent short periods, either H2 or H3 
(but never both at once) also becomes significant. These are the potentials that cause the bounces. When Hi + H2 is 
dominant at p, the evolution is essentially as if p were in one of the vacuum Bianchi II (Taub ||l8|] ) spacetimes. When 
Hi + H3 is dominant, the evolution is essentially as if p were in one of the Bianchi I with magnetic field (Rosen [p^) 
spacetimes. Both the Taub and the Rosen solutions approach one Kasner solution toward the singularity and another 
Kasner solution in the opposite time direction (r —^ —00). Given a Kasner solution there is no more than one Taub 
or Rosen solution that approaches it as r ^ — cx). So given a particular Kasner epoch at point p, one knows which 
Taub or Rosen solution approximates the next bounce toward the singularity. This allows one to approximate the 
sequence of Kasner epochs that will occur in a given Mixmaster evolution. 

To understand qualitatively why the bounces occur, let us assume that the functions P, Q, A, ttp, ttq, tt\ and their 
derivatives develop in time in such a way that they do not counteract any explicit exponential decay in any of the 
terms in H or the resulting evolution equations. We shall call this "assumption A." For example, if at some spacetime 
point, {p,t), one has the following (which we shall call "the Kasner conditions"), 
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r>0, A + T<0, P-r<0, 
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then assumption A implies that at that point Hi ^ H2 and Hi ^ H3; and further it implies that the terms in 
the evolution equations which are derived from Hi dominate those derived from H2 and H3. The relative values of 
Hi , H2 and H3 accurately monitor the relative importance of the terms derived from them in the evolution towards 
the singularity because the exponential factors control the growth and decay of the terms in which they are present. 
Without exception, our numerical results support assumption A, and we assume it throughout the following. 

Let us say that the Kasner conditions are satisfied at some point (p, t), so the evolution at (p, t) is dominated by 
Hi. If H2 and H^ were zero, the evolution at p would be exactly Kasner. The quantity which determines which type 
of bounce will occur next is v = \J Hi/ttx. One calculates (using assumption A along with the Kasner conditions) 
that dHi/dr « and dnx/dr « 0, so during a given Kasner epoch at p, v is essentially constant in time. We now 
argue that if u < 1 there will be a magnetic bounce in finite time, while if d > 1 there will be a curvature bounce in 
finite time. 

The reason v is so important during a Kasner epoch in determining the next bounce is because H3 /Hi is controlled 
by e^^'^ and H2/H1 is controlled by e^~'^, and in turn A + r and P — t are controlled by v. In particular, we find 
that during a Kasner epoch, dX/dr « — f^, so that d(A + r)/(ir « (1 — w^). Thus if u < 1, one has e^^'^ increasing 
exponentially with t and therefore H3/H1 — > 0(1) in finite time; while if w > 1, it follows that A + r decreases with r 
and H3/H1 stays small. The evolution for P, governed by dP/dr — ttp/27Tx, is a bit more complicated. However, one 
finds —V < dP/dr < +v, so if w < 1, then d{P ~ T)/dT < and H2/H1 must decrease. Hence if v < 1, a magnetic 
bounce must occur. Now consider the case t;>l. If7rp>0 and therefore dP/dr > 0, then the Kasner evolution leads 
to 7Tp/2tt\ — > V, and hence d{P — T)/dT > 0. It follows that H2/H1 — > 0(1) in finite time and there is a curvature 
bounce. If ttp < and so dP/dr < 0, then the Kasner evolution for P has a single minimum after which ttp > and 
the evolution proceeds to a curvature bounce as just noted. The evolution of P past its minimum is called a kinetic 
bounce. The kinetic bounce, caused by e~'^^7rQ/ATT\ in Hi, keeps H2 from dying off in these spacetimes. (See Figs. 

1-3.) 

What happens after a given bounce occurs? Since, as noted above, a magnetic bounce is essentially a Rosen solution 
and a curvature bounce is closely approximated by a Taub solution, one can use the known features of those solutions 
to determine the following 1) After a magnetic bounce, induced by w < 1 at p, the metric evolution at p returns to 
a Kasner epoch, this time with v > 1. A curvature bounce will eventually follow. 2) After a curvature bounce induced 
by 1 < u < 3, one returns to a Kasner epoch with t; < 1, so a magnetic bounce will follow. 3) After a curvature 
bounce induced by u > 3, a Kasner epoch occurs with w > 1, so another curvature bounce will follow. 

The preceding characterizes the behavior at one point in space. Since v, P, and A are functions of 9, nearby points 
will in general not reach the end of a Kasner epoch at exactly the same time. (See Fig. 4.) 

There are two different kinds of exceptions to the behavior just described. Neither violates assumption A, but the 
argument that H2 and H3 continue to decay and then grow again breaks down in the following ways. First, there 
exist non-generic spacetimes in the class we are considering which are not Mixmaster. For instance, if we set B = 0, 
then H3, which causes the magnetic bounces, is missing and these spacetimes will be AVTD. 

The second type of exception happens in a generic spacetime, but only at isolated points, not on an open set in 
the spacetime. There are a number of field configurations that prevent bounces at a point. Some of these have been 
seen in studies of the Gowdy spacetimes p^ : Q' = when a curvature bounce would normally occur prevents the 
curvature bounce, and ttq — when a kinetic bounce would normally occur prevents the kinetic bounce (and hence 
the next curvature bounce if w > 1). Others are new: v = 1 during a Kasner epoch prevents both the curvature and 
magnetic type of bounce and the Kasner epoch persists. If f = 3 during a Kasner epoch a curvature bounce will 
occur, but the subsequent Kasner epoch has v = 1. If vrp = during a magnetic bounce the next curvature bounce 
is prevented. While more and more of these exceptional points occur in a given spacetime as the evolution continues, 
they will, for generic initial data, always be at isolated values of 9 [ po[ . 

Our numerical study, combined with qualitative analysis of the evolution equations, provides strong evidence that, 
generically, spacetimes in this class exhibit Mixmaster behavior. While this class is spatially inhomogeneous, it is still 
very restricted. We predict that the particular topology (which determines the boundary conditions on functions of 
9) chosen is not necessary and that the generic cosmological spacetime with local 2-torus symmetry and a magnetic 
field perpendicular to the symmetry directions will also have Mixmaster behavior. But this is still a very restricted 
class. The question remains whether cosmological spacetimes in general do indeed exhibit Mixmaster behavior, which 
would be a surprising simplification of their evolution in the neighborhood of the initial singularity, or whether some 
other possibility in their evolution toward the initial singularity exists Q . 
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FIG. 1. Typical evolution at one spatial point. The thick dashed line is log]^o('if') ^^'^ the thick solid line is log]^Q(^) (left 
axis scale). The thin solid line with circles is P (right axis scale). The bounces are labeled. Note that H2 ~ Hi during the 
curvature bounces and H3 ~ Hi during the magnetic bounce. starts to grow after the kinetic bounce. 



4 



n . 6 




FIG. 2. Evolution at same spatial point as in Fig. 1. The solid line is P and the dotted line is tvq (left axis scale). The 
dash-single-dotted line is Q and the dash- triple-dotted line is ttp (right axis scale). Note that P is essentially linear in t except 
during bounces. The other functions are essentially constant for most of the evolution. This graph shows to which type of 
bounce each function responds. 
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FIG. 3. Evolution at the same spatial point as in Figs. 1 and 2. The solid line is A (left axis scale). The dotted line is tta 
(right axis scale). Note that A is essentially linear in r except during bounces. 




FIG. 4. This shows during the first curvature bounce of Figs. 1-3 with neighboring d values included, is essentially 
constant in time during a Kasner epoch and changes during a curvature or a magnetic bounce. 
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